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Abstract 

The gluon dominance model presents a description of multiparticle production in 
proton-proton collisions and proton-antiproton annihilation. The collective behav¬ 
ior of secondary particles in pp-interactions at 70 GeV/c and higher is studied in 
the project ’’Thermalization”. The obtained neutral and charged multiplicity 
distribution parameters explain some RHIC-data. The gluon dominance model is 
modified by the inclusion of intermediate quark topology for the multiplicity dis¬ 
tribution description in the pure pp-annihilation at few tens GeV/c and explains 
behavior of the second correlative moment. This article proposes a mechanism of 
the soft photon production as a sign of hadronization. Excess of soft photons allows 
one to estimate the emission region size. 


1 Introduction 

A new model of investigating multiparticle production (MR) at high energy is proposed. 
It is based on multiplicity distribution (MD) description of different interactions on basis 
of QGD and a phenomenological hadronization scheme. It is shown that the proposed 
model agrees with experimental data in a wide energy region and, perhaps, can be used 
for analysis of jet quenching and other phenomena at RHIG 

Application of this model approach to pp-interaction (for the beginning see j2]) is given 
in Section 2. The additional investigations of MD in the pp annihilation channel at a few 
tens GeV/c are carried out in Section 3. The emission region size for soft photons and 
the possible mechanisms of their formation are discussed in Section 4. The main results 
of these studies are given in Section 5. 

2 MD in interact ions (continuation) 

MD of charged particles in proton interactions by means of the gluon dominance model 
were studied in j2]. It is interesting to get MD for neutral mesons. For this purpose 
we take experimental mean multiplicity of vr® in pp-interactions at 69 GeV/c (-y/i — 
11.6 GeV). It was be found 2.57 ± 0.13 jSj. So the mean multiplicity in this process is 
calculated as the product of the mean number of evaporated active gluons (m = 2.48) and 
hadron parameter n^. We can determine the hadronization parameter for neutral mesons: 

^Talk given at XVII International Baldin - Seminar ’’Relativistic Nuclear Physics and Quantum Chro¬ 
modynamics”. JINR. September 27 - October 2, 2004, Dubna, Russia. 



n'j=1.036± 0.041 in We expect approximate equality of probabilities of different hadron 
production at the second (hadronization) stage. MD for neutral mesons have a form as 
for charged particles [2]: 
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and can be easily obtained if they are normalized to mean multiplicity tt^’s (Fig. 1). From 
this distribution we see that the maximal possible number of from TSTM jl] is 16. 
MD for the total multiplicity are shown in Fig. 2. The maximal total number of particles 
in this case is equal to 42. 

The dependence of the mean multiplicity of neutral mesons uq versus the number of 
charged particles rich can be determined by means of MD Pn^^^[s): 
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where ni and n 2 are lower and top boundaries for the total multiplicity at the given 
number of charged particles rich- The MD of charged and neutral secondaries obtained by 
TSTM give the maximal number for charged rich = 26, neutral uq = 16 and total ritot = 42. 
That is why we have the following limits for rii and 77 - 2 : rii > rich, 772 < 16 + rich- These 
restrictions result in great disagreement with experimental data [S] at small multiplicities. 
It was shown in jl]. 

A signihcant improvement will be reached if we decrease the top limit at low multi¬ 
plicities [rich < 10) to 772 = 277 c/i. This corresponds to the case when the maximal number 
of neutrals is equal to the number of charged particles, and a double excess of neutral 
mesons over positive (negative) pions is possible. Fig. 3 shows that multiplicity of neu¬ 
trals versus rich when 772 is taken equal to 2nch at small rich and 112 = IQ + rich at rich > 10. 
This restriction in m indicates that AntiCentauro events (a large number of neutrals and 
very few charged particles) must be absent. Centauro events (a large number of charged 
particles and practically no accompanying neutrals) may be realized only in the region of 
high multiplicity. 

It is assumed that at the second stage different kinds of quark pairs from the gluon 
(maximal possible number is equal to Wot) occur with equal probabilities. We will try 
to consider the formation of neutral and charged mesons as an example of the above 
assumption. The uu and dd quark pairs may appear at sufficient energy. At the end 
of hadronization the formation of two charged mesons (the law of charge conservation 
of quarks) may take place. Production of an additional neutral particle is not necessary 
while formation of a neutral meson. So we can claim that the number of charged hadrons 
will be larger than the number of the neutral ones, or the probability of the charged 
hadron production is higher than of the neutral ones. We can estimate these probabilities 
in GDM. 

MD of from one gluon at the second stage may be described by the binomial 
distribution = Cn°pQ°Pc^~‘^°. Here rit is the total number of hadrons formed from 
gluon, rif) - the number of neutral mesons among these secondaries (the number of charged 
hadrons ric = rit — tiq), Pc{po) -fhe probability of production of charged pair (one 7r°). 




The normalized condition is Po + Pc = 1- From TSTM we have obtained rich = 1-63 
and no = 1.036. The mean multiplicities for binomial distributions will be equal to: 
^ch = Pc^t, nQ = po^t- The probability of the charge particle production is higher than of 
the neutral mesons {rich > no). The ratio of these values is Pc/po ~ 1.46. 

The mean multiplicity of newly born hadrons (charged or neutral) in proton interac¬ 
tions in GDM is equal to the product of the mean multiplicity of gluons obtained at the 
first stage and the mean multiplicity of hadrons {rich Fq) produced from one gluon at 
the second stage. In the case of binomial distribution rich = nfPc, no = rifPo. Taking into 
account two leading protons, the mean multiplicity is nch{s) = 2 + rrig{s) - rich for charged 
particles in pp-interactions. The mean multiplicity of neutral mesons in this process is 
no(s) = nig{s) ■ Hq. The ratio of the mean charged pairs to the neutral mesons in proton 
interactions is 

nch{s)/‘2 ^ 1 ^ 1 . ^ /gN 

no m(s) • nQ 2 

At 69 GeV/c this ratio Q is equal to 1.19 ± .25. At the higher energy the mean number 
of active gluons m increases and becomes much more than 3. In this case 0 it will be 
around the ratio of n^i^/2nQ. The experimental data have shown 1.6 for Au-Au peripheral 
interactions (80 — 92(%) centrality class) at 200 GeV and for pp interactions at 53 GeV 
jnj. We can compare these results with GDM at higher energies. 

The application of GDM to describe MD in the energy region (102, 205, 300, 405 and 
800 GeV/c) j7] in both schemes (TSMB and TSTM) j2] leads to good results (Fig. 4-8). 
Parameters of TSTM in this domain are given in Table 2. 

Table 1. Parameters of TSTM. 
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We see that the number of active gluons and their mean multiplicity increase, pa¬ 
rameters of hadronization N and n\ vary very slowly. At these energies the charged 
hadron/pion ratio Q grows up to 1.6. The parameter of hadronization has a trend 
to increase weakly but rig does not almost change. This behavior may be related with the 
production of other charged particles (not only pions): protons, antiprotons, kaons and 
so on. We consider that parameter ri^y^ goes to the limit value (like saturation). 

On the other side a small growth rig in proton interactions also points at a possible 
change mechanism of hadronization of gluons in comparison with the transition gluons 
to hadrons in e’''e“ annihilation. It is considered that in the last case partons transform 
to hadrons by the fragmentation mechanism at the absence of the thermal medium. Our 
MD analysis gives ~ 1 for this fragmentation [Hj. The recombination is specific for the 
hadron and nucleus processes. In this situation a lot of quark pairs from gluons appear 
almost simultaneously and recombine to various hadrons [H]. The value rig becomes 
bigger ~ 2 — 3), that indicates to the transition from the fragmentation mechanism to 






















the recombination one.The recombination mechanism provides justihcation for applying 
the statistical model to describe ratios of hadron yields (the ratio Baryon/Meson ~ 1) 
inj. The collective flow of quarks may be explained by the recombination mechanism, 
too. The rapid local thermalization may be a consequence of this formation of secondary 
hadrons p. 

In this way we try to compare two kinds of processes which have different values of 
hadronization parameters. The hrst one is e''“e“ - annihilation. It is usually supposed that 
fragmentation dominates in it and newly formed hadrons fragment with a high moment 
of parton into the surrounding vacuum (such objects can also appear from the hot surface 
in peripheral events in nucleus and hadron collisions) p. 

The nuclear modihcation factor Rcp and elliptic flow V 2 in Au-Au collisions at RHIC 
have revealed an apparent quark-number dependence in the pt region from 1.5 to 5 
GeV/c. Moreover, the baryon production increases more rapidly with centrality than the 
meson production. These observations conhrm the picture of hadron formation by quark 
recombination P and point out that the hadronization processes in high energy nucleus 
interactions are modihed to the comparison of e"*" -1- e~ and partly p + p collisions. 

The GDM with a branch gives growth of the part of the evaporated gluons to 0.85-0.98 
and a small rise of gluon branch number at higher energies. Besides we have got data 
about emergence of hard constituent in MP QDI. In GDM it can be explained not only 
by not only evaporation of a single gluon sources but also of groups with several gluons 
(formed by branch). A simple MD scheme of this superposition will be analyzed below. 

Let us compare MD (P with the descriptions of experimental data obtained by various 
approaches. We bring two of them. A fortunate expression for KNO function was obtained 
by a group from IHEP m who combined the elastic and inelastic processes. We can see 
(Fig. 9) good agreement with data P at 800 GeV/c both of MD in MGD (solid line) and 
KN0-function (dot line). 

A wide research of MD in pp-interactions was fulhlled by L.Van Hove, A.Giovannini 
and R.Ugoccioni They proposed a two-step mechanisms of MP. The independent 
(Poisson) production of groups of ancestor particles (named ’’clan ancestors”) were supple¬ 
mented by their decay, according to a hadron shower process (the logarithmic MD within 
each clan). Such convolution of two mechanisms gives a negative binomial distribution 
(NBD) for hadrons 


^ ^ ^ _ h{kh + l) ...jh + n-l) f n{s) Yf kh .. 
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where kh - the NBD parameter and n{s) - the mean multiplicity of hadrons. The compar¬ 
ison of NBD (dot line) and our MD in GDM (solid line) with data at 800 Gev/c is given 
in Fig. 10. A. Giovannini emphasizes that the nature of this clan is gluon bremsstrahlung 
HU- Our investigations by GDM allows to give a concrete gluon content. Binomial distri¬ 
butions (BD) describe the hadronization stage. The clan model of ^U] uses the logarithmic 
distribution of secondaries in a single clan. Both of MD have the similar behavior. 

At the top energy (especially at 900 GeV) the shoulder structure appears in [T^ . 
The comparison of data with one NBD does not describe data well. But the weighted 
superposition of two NBD gives a good description of the shoulder structure is) HD]. 
At 14 TeV A.Giovannini expects the weighted superposition of the three classes of events. 






Figure 1: MD for 7r°. 


Figure 2: MD for ritot- 


Figure 3: nlijich)- 


We can modify our GDM considering that the gluon fission may be realized at higher 
energies. The independent evaporation of gluons sources of hadrons may be realized by 
single gluons and also groups from two and more fission gluons. Following A.Giovannini 
we name such groups of gluons - clans. Their independent emergence and following 
hadronization content of GDM. MD in GDM with two kinds of clans are: 
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where ai and a 2 are the contribution single and double gluon clans (ai + 0:2 = 1). The 
comparison with experimental data for proton interactions at y/s = 62.2 GeV na is 
given in Fig. 8. We have obtained the following values of parameters: N = 7.06 ± 3.48, 
mi = 3.59 ± 0.03, m 2 = 1.15 ± 0.25, Uh = 3.23 ± 0.14, Mgi = 8, Mc /2 = 4, ai/aa ~ 1.8 
at x^/ndf=9.12/13. The mean multiplicities of the two kinds of clans are similar. 

The specihc feature of our GDM approach is the dominance of a lot of active gluons 
in MP. We can expect the emergence of them in nucleus collisions (experiments at RHIG) 
and the formation of a new kind of matter (quark-gluon plasma)at high energy. We 
consider that our gluon system can be a candidate for this. So the mean multiplicity of 
active gluons approached 10 at RHIG. For Au-|-Au central collisions their number may be 
equal to 200 ■ m 2000 before the branch. This gluon medium facilitates the quenching. 


3 MD in 2^]9-annihilation 

In the midst of interesting and enough inextricable hadron interactions the pp annihilation 
shows up especially HH. Experimental data at tens GeV/c ^3] point out on some maxima 
in differences between pp and pp inelastic topological cross sections what may witness 
about the contribution of different mechanisms of MP 


Act„(PP - pp) = <Tu{PP) - ffn(PP)' 


( 6 ) 











The important information about the MP mechanism may be picked out from the MD mo¬ 
ment analysis of charged particles. The second correlative moment for negative particles 
/2 are available to study MP 

f 2 ~ =n-{n- -1) -nJ. (7) 

The negative value of second correlative moments is characteristic for a more narrow 
MD than Poisson, and they indicate the predominance of the hadronization stage in MP. 
According to MGD, active gluons are a basic source of secondary hadrons. 

At the initial stage of annihilation three valent gg-pairs {uud and uud) are. They 
can turn to the ’’leading” mesons which consist from: a) valent quarks or b) valent and 
vacuum quarks [S] . In the case a) only three ” leading” neutral pions (the ” 0” topology) or 
two charged and one neutral ” leading” mesons (” 2” - topology) may form. In b) case the 
”4”- and ”6”- topology is realized for ’’leading” mesons. We suggest that the formation 
neutron and antineutron (exchange) can be realized. 

A simple scheme of MP for annihilation may give the negative second correlative mo¬ 
ments in GDM. We suggest that the active gluon emergence together with the formation of 
intermediate topology occurs. The GF for a single active gluon Qi(z) = [l-|-n/iV( 2 ; —1)]-^ 
gives jlj 

/2 = Ql(z)U=i - [Qi(;^)|.=i]' = -{n^f/N < 0. (8) 

Reciprocally for m gluons GF and /2 will be 

(3„(z) = [l+il‘/iV(z-l)]"", U = -m(rt‘flN. (9) 

We consider that m grows while increasing the energy of the colliding particles, and /2 will 
decrease almost linearly from m. Such behavior qualitatively agrees with experimental 
data 113]. If we take concrete MD for gluons, then GF for secondary hadrons and /2 
are 

Q(.A ^iN(z - i)r" (10) 

m 

f 2 = [f 2 = m{m-l) (11) 

where - the second correlative moment for gluons. In this scheme /2 may be negative or 
positive. We consider that the negative value /2 in the large energy region in comparison 
with p + p interactions may be related with the destruction of the initial system on three 
or more shares and the number of active gluons related with a ’’leading” pion will be less 
than in the case of a leading proton in pp-collisions at the same energy. Herewith the 
total number of such gluons at annihilation may be bigger, their manifestation happens 
independently but the number of them per one pion grows slowly. The explanation of 
the negative /2 was given R.Lednicky usi at the assumption of the independent MP of 
charged particles. The second correlative moment has a zero value only in the small 
energy domain. And so we should restrict the region to apply this explanation. 

According to GDM for pp annihilation and taking into account three intermediate 
charged topology and active gluons, GF Q{z) for final MD may be written as the convo¬ 
lution gluon and hadron components: 







Figure 4: MD at 102 Figure 5: MD at 205 Figure 6: MD at 303 
GeV/c. GeV/c. GeV/c. 
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Figure 7: MD at 405 Figure 8: MD at 800 Figure 9: MD in GDM and 

GeV/c. GeV/c. KN0-function [TT] . 


The parameters of cq, C2 and C4 are determined as the part of intermediate topology (”0”, 
”2” or ”4”) to the annihilation cross section (cq -|- C2 -t- C4 = 1). For the simplicity we are 
limited by Poisson distribution with the finite number of gluons for P^. 

The comparison of the experimental data (Fig. 12) gives the following values of pa¬ 
rameters: m = 3.36 ± 0.18, N = 4.01 ± 0.61, n^=1.74 ±0.26, the ratio cq : C 2 : C 4 = 15 
: 40 : 0.05 at = 5.77/4 and the maximum possible number of gluons M = 4 at 

”4”-topology. The sum begins from m = 1 (inelastic events), at n > 2 - from m = 0 and 
finishes up m < M at small multiplicities (n < 4). We should to emphasize very compli¬ 
cated events rich = 0 and 2. This research of pp annihilation requares to be continued. 
We will develop MGD to describe MD at energies 200, 500, 900 GeV jlB] and higher. 


4 Soft photons 

The production of photons in particle collisions at high energies was studied in many 
experiments irz|. In project ” Thermalization” it is planned to investigate low energetic 
photons with pt < O.lGeD/c and x < 0.01 [TH]. Usually these photons are named soft 
photons (SP). Experiments shown that measured cross sections of SP are several times 
larger than the expected ones from QED inner bremsstrahlung. Phenomenological models 
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Figure 10: MD in TSTM 
and NBD. 


Figure 11: MD in GDM 
(clan). 


Figure 12: GDM MD in pp 
at 14.75 GeV/c. 


were proposed to explain the SP excess: the glob model of Lichard and Van Hove and the 
modihed soft annihilation model of Lichard and Thomson ra • 

We consider that at a certain moment QGS or excited new hadrons may set in an 
almost equilibrium state during a short period or hnite time. That is why, to describe 
massless photons, we will try to use the black body emission spectrum m- From ex¬ 
perimental data [IH] the inelastic cross section is equal to approximately 40m&, the cross 
section of SP formation is about Arab, and since ~ ■ ain, then the number of SP 

will be equal to ~ 0.1. For convenience, we may use the well-known density of MVB at 

Tr = 2.27577 and get the number of photons by means of MVB n-y{T) = n^{Tr) ■ . 

The density of SP in the region 1/m^ will be equal to 

p{T) = n^{T)/V = 4.112 ■ 10® • 10“® • 10“®^ ■ fmr^. 

The estimates of temperature are implemented by transfer moment: T = p ^ PtV^ 
{IMeV = 1.16 ■ 10®°77). If T{pt) is known, using we can estimate the linear size of 
radiation system {V ~ L®). Dependencies of the linear size of system (L) from the SP 
moment {pt) are given in Table 2. 

Table 2. The size of system L (fm) versus pt (MeV/c) of SP. 


Pt 

10 

15 

25 

30 

40 

50 

L 

11 

6.9 

4.1 

3.5 

2.6 

2.0 


It is well-known that the temperature of second hadrons is higher than the temperature 
of SP. We presume that objects with soft gluon content may not transform into hadrons 
but turn into SP. The amount of such soft gluons is estimated by Ng in TSMB. 


5 Conclusion 

In our research we have undertaken an attempt to give MP description in different pro¬ 
cesses by means of a unihed approach based on quark-gluon picture using the phenomeno- 
















logical liadronization model. The implemented model investigation allows us to under¬ 
stand deeper the picture of MP at various stages. We have obtained qualitative and 
quantitative agreements of our schemes with experimental data in e+e”, pp annihilation 
and pp and nucleus collisions in a very wide energy domain. 

The authors appreciate for the support of physicists from JINR, GSTU who encour¬ 
aged our investigations. Project ”Thermalization” is partially supported by RFBR grant 
03-02-16869. 
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